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Abstract 

We investigate the efficiency of single timeslice stochastic sources for the calculation of 
light meson masses on the lattice as one varies the quark mass. Simulations are carried out 
with Nf = 2 flavours of non-perturbatively 0{a) improved Wilson fermions for pion masses in 
the range of 450 — 760 MeV. Results for pseudoscalar and vector meson two-point correlation 
functions computed using stochastic as well as point sources are presented and compared. 
At fixed computational cost the stochastic approach reduces the variance considerably in the 
pseudoscalar channel for all simulated quark masses. The vector channel is more affected by 
the intrinsic stochastic noise. In order to obtain stable estimates of the statistical errors and 
a more pronounced plateau for the effective vector meson mass, a relatively large number of 
stochastic sources must be used. 
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1 Introduction 



In lattice QCD, hadronic properties such as masses and matrix elements can be computed in 
terms of Euclidean correlation functions. Typically, these are expectation values of properly 
selected polynomials in the quark and gluon fields that project on states with the desired 
quantum numbers. After performing the Wick contractions, correlation functions are ex- 
pressed as traces over products of quark propagators, Dirac matrices and color-structures. In 
order to obtain precise estimates of hadron masses and transition amplitudes, variance reduc- 
tion methods may be applied such that correlation functions with good statistical accuracy 
can be computed. 

The quark propagator in coordinate space can be computed as the solution of the linear 
system 

D^ = r,, (1) 

where D is the lattice Dirac operator and rj a source vector. In its simplest form, rj is taken 
to be a point source, i.eQ 

7?(X') = 6x>y . (2) 

This implies that the solution of eq. ([T]) yields the quark propagator from a single point y to 
any other point x, which corresponds to just one column of the propagator matrix. Thanks 
to translational invariance a large class of correlation functions can be defined in terms of 
these "one-to-all" propagators. However, in this way only a small fraction of the information 
contained in D is processed. In typical simulations of lattice QCD the sparse matrix D 
has 0(10^ X 10^) entries, and solving eq. ([!]) to machine precision for all source positions is 
therefore beyond the capabilities of even the most powerful supercomputers. Volume-filling 
random-noise sources have been proposed as a means to access the full propagator matrix [l}|4] 
by replacing it with a stochastic estimate. These stochastic "all-to-all" propagators have been 
successfully applied in a number of different contexts (2| jTH[T4|[l5] . However, care must be 
exercised in their actual construction, since an arbitrarily chosen stochastic source vector can 
induce a large variance into hadronic correlation functions. In other words, for the method 
to be efficient, one must ensure that the intrinsic stochastic noise does not overwhelm the 
gain in information provided by having access to the entire propagator matrix. A particular 
technique, which proved to be quite efficient, is the so-called "one-end trick", pioneered in [Sp] 



and successfully applied in several studies of light meson physics 10, 11, 14, 15 . 

In this article we present a systematic study of the effectiveness of single-timeslice stochas- 
tic sources. In particular, we monitor the variance and the computational cost as a function 
of the quark mass and compare it to the variance achieved with point sources. Our simu- 
lations are performed for QCD with A'^f = 2 flavours of 0{a) improved Wilson quarks. We 
concentrate on two-point correlation functions of the pion and the rho meson. 

The outline of this paper is as follows: in section[2] we revisit the concepts of stochastic 
volume sources and the one-end trick. The set-up of our simulations is discussed in section|3j 
and our main results are presented in section|4j Finally, section[5] contains a summary of our 
findings as well as some concluding remarks. 



^For simplicity we suppress colour and spinor indices. 
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2 Stochastic propagator estimation 



In this paper we restrict ourselves to two-point correlation functions of a flavour off-diagonal 
quark bilinear, Oudi^) = 'u{x)Td{x), where u and d denote the fields of the up- and the down- 
quarks, respectively, and T is one of the 16 Dirac matrices specifying the desired quantum 
number. The two-point correlation function is given by 

(^OUx)OUy)) = (Tv[s{x,y)TS{y,x)t]y (3) 

where T = 7or"l'7o, and (...) denotes the average over the gauge configurations. Since we 
assume exact isospin symmetry, the same symbol, S, is used to denote the propagator for 
both quark flavours. The spectrum of the particles with the prescribed quantum numbers 
can then be determined from the exponential decay of the zero-momentum projection, i.e. 

C2(t;y) = ^(Tr{5(x,y)r5(y,x)f}), t = - yo- (4) 

X 

For a given gauge field the propagator S{x,y) from lattice sites y to x, is obtained as the 
solution of the linear systerrj^ 

^ ^ Z)- (x, z)S!^'^{z, y) = S^^6''%y . (5) 

Using point sources amounts to solving the linear system for one particular choice of y. 
This implies that altogether 12 inversions must be performed, since four spinor and three 
colour components must be considered independently. While this is sufficient to compute the 
correlation function C2{t; y), the statistical signal could be further improved by averaging over 
many different source points. An improved estimator for the correlation function is obtained 

by 

y x,y 

where V3 denotes the number of lattice sites within one timeslice. In order to see how such 
an average can be effected, we revisit the method of stochastic noise sources in the following. 

2.1 General stochastic noise sources 

In the stochastic approach, an ensemble of N-i- random vectors, {r/^^'-'(xo, = 1, . . . , A^^^}; 
is generated for each gauge configuration. The number Ar is sometimes referred to as the 
number of "hits" . These source vectors are created by assigning independent random numbers 
to all components, i.e. to all lattice sites, colour and Dirac indices. Each random number is 
drawn from a distribution T) which is symmetric about zero in the hit limit A^ — )• 00, i.e. 

(C(^o,f)>src ^ lim - (^M):(xo,f) = 0. (7) 

r=l 

In addition the sources satisfy 

(<(f,Xo)(r/"f)^(y,2/o)) = ^xoyo^xy^a^S''^ ■ (8) 
\ I src 

^ Latin indices are used to label color, while Greek letters denote spinor components. 
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Solving the linear system of eq. ([T]) for each of the N,- source vectors yields a set of solution 
vectors 

(ci>w);^(x) = Y.H ^) iv^^i^)- (9) 

z c,7 

The estimate of the entire propagator is defined as the stochastic average ( "hit average" ) over 
the product between solution and noise vectors 



z c,7 

t remains to specify the distribution from which the random vectors 77^'"^ are drawn. In refs. 



l][2 it was noted that a flat distribution of Z(2)-elements, (r?^'"))!^ (x) e V = Z(2) = {+1, -1}, 
or, more generally, elements oiT) = 'L{N) is very effective in realising the condition of eq. ([8|. 
In this work we follow Foster and Michael [S] and draw separate elements of Z(2) for the real 
and imaginary parts of the source vector, i.e. 

P = Z(2)$DZ(2) = |^(±l±i)|. (11) 

Experience shows that a random source vector, which is distributed over the entire space-time 
lattice, leads mostly to a very noisy signal for hadronic correlation functions. An essential 
step towards a significant variance reduction is taken by restricting the support of the source 
vector to individual timeslices, Dirac or colour components js]. Such "dilution schemes", 
and, in particular, the so-called time-dilution are widely used in the computation of hadronic 
properties [5|[7||8]. Here, the non-zero components of the random source vector are restricted 
to a single timeslice, e.g. yo = 

"^^y' " \ 0, otherwise ' ^ > 

We end this discussion with a few general remarks. In practice the limit A^r — )• 00 cannot be 
taken, and thus the details of constructing the random source should be optimized for the 
correlation function at hand. In ref. |11| it was noted that for many correlators the stochastic 
and gauge averages commute. It then suffices to generate a reasonably small number of 
random source vectors per gauge configuration, since the ensemble average automatically 
implies the stochastic average. 

2.2 The one-end trick revisited 

The naive implementation of stochastic sources consists in replacing each propagator in the 
correlation function of eq. ([g]) by the stochastic estimate of eq. (10). For a general choice of 



F this implies that an independent source vector must be used for each quark propagator, 
in particular if F couples different Dirac components. This typically results in a relatively 
noisy signal [12] . By contrast, for correlators which involve only diagonal combinations of 
spinor components, it is possible to compute the two-point function stochastically using only 
a single random source, which is distributed over all colour and Dirac components within one 
timeslice. The relative numerical effort compared to using a point source is thereby reduced 
by a factor 12 per hit: Setting Nj- = 12 results in the same number of inversions that must 
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be performed, while the statistical error can be expected to be reduced. This is the so-called 
"one-end trick" (5[[9]. 

The situation in the case of correlation functions, in which different spinor components are 
coupled, is less favourable but can be dealt with via the generalized one-end trick, sometimes 
also referred to as the "linked source method" |14|. It amounts to choosing a spin-diagonal 
random source vector, which has support only for a particular spin component r and a single 
timeslice yo (e.g. yo = 0), viz. 

{r]^'Kiy) = {C^'^)\y)SoyoS.r, (13) 

where the components of the stochastic vector ^^''^ are drawn from a distribution T> and satisfy 

(e{m^)\y)) =6,yS^'. (14) 

\ / src 

Solving the linear system of eq. ([T]) for spin component r yields the solution vector ^{x), i.e. 

y b 

The correlation function at zero momentum is then obtained as 

^2(t) = -(EE (f(r75)'&(:r)tj" [(75f)cl>(x)j" \ ), (16) 

which is the stochastic estimator of the two-point function at zero momentum of eq. ^ after 
applying eq. Q. The generalized one-end trick can be applied whenever the stochastic source 
vector commutes with the given choice of F-matrices. The spin-diagonal source vector ^^""^ 
satisfies this requirement by construction. Compared with the point source, the numerical 
effort is reduced by a factor three per hit. 

The use of linked sources is not mandatory for pseudoscalar mesons, since the diagonal 
F-structure of the associated correlators (F75 = 1) implies that stochastic noise vectors com- 
mute without any modification. However, during the course of a simulation, many different 
hadronic channels, involving arbitrary Dirac structures at both the source and sink, are con- 
sidered. In order to facilitate the cost comparison for different correlator channels, we have 
implemented linked sources by default, and below we proceed to compare their effectiveness 
in simulations covering a range of light quark masses for both the pseudoscalar and vector 
channel. According to ref. fTl|, the use of linked sources for pseudoscalar mesons was not 
found to be inferior compared with non spin-diagonal sources, at fixed computational cost. 

One potential drawback of the one-end trick is that two-point functions cannot be com- 
puted for arbitrary momenta using a given set of noise vectors. Owing to the automatic 
summation over the spatial source and sink coordinates, a specific momentum is selected. To 
utilize the one-end trick for computing the two-point correlator at a given non-zero momen- 
tum, the set of noise vectors {i]^^\yo,y)\r = 1, . . . ,NrJ must be transformed separately for 
each selected momentum with an appropriate phase e*^^ prior to performing the inversion. It 
is worth noting that (partially) twisted boundary conditions |13j can be successfully combined 
with random sources [14[|16] . Again, for each momentum channel, N,- extra inversions are 
then required, and thus the numerical cost increases in relation to the point source, where 
the two-point function can be projected at least on the Fourier modes at negligible additional 
cost. 
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3 Simulation setup 



This work is based on gauge configurations with N{ = 2 flavours of non-perturbatively 0(a)- 
improved Wilson fermions which have been generated as part of the CLS effort |17j, us- 



ing the deflation-accelerated DD-HMC algorithm 18 . This algorithm combines domain- 
decomposition (DD) methods p.9| with the Hybrid Monte Carlo (HMC) algorithm (20| and 
the Sexton- Weingarten multiple-time integration scheme |21j . All ensembles in this work were 
generated for f3 = 5.3, a choice for which the coefficient of the Sheikholeslami-Wohlert term 
was determined as Csw = 1.90952 [22]. 



Run 


Lattice 


Number cfgs. 


^sca 




[MeV] 


E2 


64 X 323 


158 


0.13590 


0.24292(29) 


760 


E3 


64 X 32^ 


156 


0.13605 


0.20645(37) 


645 




64 X 32^ 


162 


0.13610 


0.19305(41) 


605 


E5 


64 X 32^ 


159 


0.13625 


0.14345(55) 


450 



Table 1: Simulation parameters and pion masses. The latter are determined using the one-end-trick 
with = 6 hits each on three different timeslices for the source vector. 

Our simulation parameters are listed in Table[T] The ensembles contain configurations 
which are sufficiently separated in Monte Carlo time such that autocorrelations can be 
safely ignored. In this work we always set the valence quark mass equal to that of the 
sea quarks. The inversions of the Wilson-Dirac operator of eq. ([T]) were performed using a 
Schwarz-preconditioned generalized conjugate residual (SAP-I-GCR) algorithm [23]. Further 
simulations details are described in [24| . The conversion of pion masses into physical units was 
performed using the preliminary scale determination via the mass of the Omega baryon ^25j , 
which yields a = 0.063 fm at /3 = 5.3. 

Mesonic two-point correlators were computed in the pion (F = 75) and rho (F = 7j) 
channels. For the latter we averaged the contributions from all three spatial 7-matrices. 
Only flavour non-singlet correlators were considered. The quark propagators entering the 
correlation functions were computed either using a point source or by applying a random 
source in the manner of the generalized one-end trick described above. For the latter linked 
sources with A'j- = 1, 3 and 6 hits were used. While for N,- = 3 hits the numerical effort 
expressed in terms of the number of inversion remains the same as for the point source, it 
reduces by a factor of three for Nj- = 1. Six hits represent twice as many inversions compared 
to the point source. In order to study the scaling of the variance with A^r in a more detailed 
manner, we used as many as 25 hits for ensemble E4. The accuracy of our calculations is 
enhanced by averaging the results obtained for three different locations for each type of source, 
corresponding to source positions Xo/a = 0, 21 and 42. 

In our analysis the forward-backward symmetry of the correlators was exploited to average 
the data about the central timeslice T/2. In order to extract the masses of the ground 
state we performed correlated x^-minimizing fits to the folded correlation functions. The 
next-to- lowest state was explicitly taken into account, by substituting its energy by 3m. 
in the pseudoscalar channel and 2y^m^ -|- (27r/L)2 in the vector channel, respectively 26 



Statistical errors were determined via the single-elimination jackknife method. As will be 
explained below, the achieved statistical precision on spectral quantities enters our definition 
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of a measure for the efficiency of stochastic sources. In order to quantify the significance of 
this measure we have used Berg's proposal of a nested jackknife procedure (also referred 
to as second-level jackknife) to estimate the fluctuations of the statistical error (i.e the "error 
of the error" ) . 

4 Results: Pseudoscalar and vector meson two-point correla- 
tors 

Our objective is to compare the performance of point sources and the generalized one-end 
trick when computing light meson masses. Therefore, for a variety of different quark masses 
the resulting variances are monitored as a function of the computational effort, which is 
expressed in terms of the number of inversions, A^invi which are required to solve the linear 
system of eq. ([T| for one particular source position. This is motivated by the observation that 
the number of iterations of the deflated SAP-I-GCR solver does not depend on the source type 
but only on the simulated quark mass. Note that A^inv represents the number of inversions 
per source position. 

Table[2] contains the results of the pion and rho meson mass fits using both source types 
for the various gauge field ensembles whose simulation parameters are listed in Table[T] By 
default, results from the three source positions are averaged over. The results are illustrated 
in Fig.[l| which shows the ratio of errors obtained using stochastic volume sources and point 
sources, respectively. 



Source type 


iVinv 


E2 






E^ 










an, 






Point 
Volume: 
Volume: 
Volume: 


1 hit 
3 hits 
6 hits 


12 
4 
12 
24 


0.24216(52) 
0.24337(35) 
0.24309(30) 
0.24292(29) 


0.20647(52) 
0.20684(45) 
0.20650(40) 
0.20645(37) 


0.19240(69) 
0.19335(47) 
0.19320(42) 
0.19305(41) 


0.14264(81) 
0.14402(59) 
0.14346(57) 
0.14345(55) 










an 






Point 
Volume: 
Volume: 
Volume: 


1 hit 
3 hits 
6 hits 


12 
4 
12 
24 


0.3782(38) 
0.3819(57) 
0.3843(39) 
0.3858(29) 


0.3501(49) 
0.3523(80) 
0.3518(52) 
0.3559(39) 


0.3353(56) 
0.3166(80) 
0.3367(54) 
0.3333(36) 


0.2842(80) 
0.2909(139) 
0.2820(73) 
0.2842(61) 



Table 2: Pion (upper half) and rho meson (lower half) masses computed using point sources and 
the generalized one-end trick with Nr = 1,3 and 6 hits for the ensembles E2 — E^. Ni^v denotes the 
number of inversions per source position. Fit ranges were chosen as 10 <xo/a< 30 and 11 < a;o/a < 30 
in the pion and rho meson channels, respectively. In case of the lightest rho meson the fit range was 
reduced to 11 < xo/a < 27, due to strong fluctuations around the central timeslice. 

Our investigation of the uncertainty in the statistical error estimate itself via a nested jack- 
knife procedure revealed that the fluctuations in the error estimate are quite small, amounting 
on average to about (4 — 7)% in the pion channel and (7 — 10)% in the rho meson channel, 
respectively. 

Comparing the conventional point source to the generalized one-end trick in the pion 
channel we observe that random noise sources lead to a considerable improvement of the 
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pion channel 
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Figure 1: The statistical error of meson masses computed using volume sources normalized by the 
error obtained with the point source, plotted against the squared pion mass for TVr = 1, 3 and 6 hits. 
Error bars (where shown) were obtained using a nested jackknife procedure. 
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hits 




^4 




hits 


iVinv 


E4 










arup 






am.,^ 


amp 


pt 


12 


0.19240(69) 


0.3353(56) 


8 


32 


0.19302(41) 


0.3370(35) 


1 


4 


0.19335(47) 


0.3166(80) 


9 


36 


0.19304(41) 


0.3361(34) 


2 


8 


0.19326(44) 


0.3352(67) 


10 


40 


0.19299(41) 


0.3363(33) 


3 


12 


0.19320(42) 


0.3367(54) 


13 


52 


0.19291(41) 


0.3377(32) 


4 


16 


0.19300(42) 


0.3288(43) 


16 


64 


0.19286(40) 


0.3399(29) 


5 


20 


0.19302(42) 


0.3328(41) 


19 


76 


0.19289(40) 


0.3408(29) 


6 


24 


0.19305(41) 


0.3333(36) 


22 


88 


0.19291(40) 


0.3417(28) 


7 


28 


0.19301(41) 


0.3346(34) 


25 


100 


0.19293(40) 


0.3430(27) 



Table 3: Variance reduction of the generalized one-end trick as a function of the number of hits. Pion 
and rho meson masses for = 1, . . . , 25 hits computed on the ensemble £'4. The values of the point 
source (labeled by 'pt') are taken from Table[2j 

statistical signal. Already a single random noise vector, corresponding to one third of the 
relative computational effort leads to a significantly reduced variance. The variance saturates 
very quickly when increasing the number of hits. This implies that the gauge noise dominates 
over the stochastic noise of random volume sources. Furthermore, a considerable improvement 
of at least 25 percent at equal cost is observed for all simulated quark masses. Our confidence 
in the results is supported by the quality of plateaus of effective masses. The plots for the pion 
are shown in Fig. [5] of Appendix[A| and illustrate that outliers of the plateau are suppressed 
due to the volume averaging effect of our stochastic sources. 

In the vector channel the generalized one-end trick is less effective. A single hit does not 
suffice to reach the precision of the point source method and no reliable estimates for the 
errors of the error were obtained. Therefore, in the right panel of Fig. [T] the errorbars of the 
results of a single hit are suppressed for small quark masses. However, at equal computational 
cost the variances are comparable and, in particular, for the two lightest quark masses the 
volume source slightly gains over the conventional point source. Contrary to the pion channel, 
it is seen that the stochastic noise is not immediately saturated and increasing the number of 
hits reduces the variance significantly]^ In order to study the scaling of the variance in more 
detail we have computed the pseudoscalar and vector correlators on the E4 ensemble for up 
to 25 hits. The results are shown in Table[3]and Fig.[2j 

The results compiled in the table demonstrate that the variance in the pion channel is 
saturated already after performing three hits, indicating that this correlator is dominated by 
the gauge noise. Increasing N-^ from 3 to 25 results in a small additional reduction of the 
error in the pion mass of only about 5%. 

By contrast, the stochastic noise dominates in the rho meson channel. Increasing 
from 1 to 7 produces a reduction of the statistical error by 60%, and a further 10% can be 
gained if is as large as 25. Assuming that the error scales hke 0{1/ y/Wr), one expects that 
the squared error times the number of inversions is constant. As the right panel in Fig. [2] 
shows quite clearly, this is indeed true in the vector channel for < 7. When larger values 
of Nj: are considered, the rate of error reduction slows down relative to the extra number of 
inversions, and thus the procedure becomes increasingly ineffective in terms of computational 



The effective mass plots of the rho meson are shown in Fig. [ejof Appendix [b] 
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rho channel 
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Figure 2: Effectiveness of stocfiastic sources as a function of the number of hits, N^-, for ensemble 
E4 in the pseudoscalar (left panel) and vector (right panel) channels. Solid black points represent the 
statistical error relative to the single-hit volume source. Red squares denote the variance scaled by 
the number of hits. 
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point 


Nr = 3 


iVr = 6 


Nr = 9 


Nr = 16 


point 





0.08 


0.21 


-0.10 


-0.30 


iVr = 3 


-0.07 





0.13 


-0.16 


-0.35 


Nr = 6 


-0.17 


-0.11 





-0.25 


-0.42 


iVr = 9 


0.11 


0.19 


0.34 





-0.23 


Nr = 16 


0.43 


0.54 


0.73 


0.29 






Table 4: The relative efficiency Cij in the vector channel for ensemble E4. 



overhead. 

In order to formulate a more quantitative criterion for the performance of point and 
stochastic sources, we define the relative efficiency, e^j, of two procedures i and j, as 

[(variance) X (number of inversions) ]j ^ 
[(variance) x (number of inversions)]^ 

Thus, the ratio of the squared error is scaled with the ratio of the computational cost. In 
Table[4]we list the values for eij in the vector channel for ensemble E4. The information in 
this matrix-like table must be interpreted in the following way: a negative entry eij in row i 
and column j means that procedure i is more efficient than j by (100 x |ejj|) %. From the 
numbers in Table[4] one concludes that, in the vector channel, a volume source with Nr = 3 
is only slightly more efficient than the point source, even though the number of inversions 
performed in both cases is the same. However, if one demands that the statistical error be at 
least as small as for the point source, a larger number of hits is more favourable: According 
to Table[3l Nr = 6 seems to be the optimal choice, since the extra numerical effort is more 
than compensated by the resulting reduction in the variance. 

In Fig. [3] the relative efficiency defined in eq. (17) with respect to the point source is 



plotted as a function of the squared pion mass. For a single hit the volume source is about 
a factor of five more efficient than the point source in the pion channel whereas at identical 
computational cost the volume source still outperforms the point source by a factor of two. 
As pointed out before, this loss of efficiency is due to the fact that in the pion channel the 
variance is completely dominated by the gauge noise. The gain in the vector channel is not 
so obvious. A stochastic volume source with A^'r = 6 hits appears to be a good compromise 
between numerical effort and statistical accuracy across the entire mass range under study, 
despite the large uncertainties in the relative efficiency Cij in the vector channel (see Figs.[l] 
and[3]). 

Another beneficial effect of using random sources in the vector channel can be seen from 
the effective mass plot shown in Fig. [4} By comparing the data obtained using Nr = 6 and 
Nr = 25 hits with those of the point source, one clearly sees that not only the statistical error 
decreases for a large hit number but that the overall quality of the plateau is much improved 
as well. This should make for much more reliable estimates of the mass in the vector channel. 

Finally we investigate the effects on the variance of distributing random noise sources 
on several timeslices instead of placing them on a single one. Let Ntot denote any number 
of inversions of the Dirac operator. If the random noise source is placed on Nts different 
timeslices, and if Nr hits are performed for each of these source positions, then Atot = Nts-Nr. 
Can the variance be reduced by choosing Ats and A^ such that Atot stays fixed? 
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Figure 3: Relative efficiency as defined in eq. (17) witli respect to tfie point source as a function of 



the squared pion mass for the pion (left panel) and rho meson (right panel) using = 1,3 and 6 hits. 
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Figure 4: Effective masses in the vector channel obtained by means of the point source and random 
sources with 6 and 25 hits. 
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iVts = 1 [0] 


Nr 


iVts = 2 [0,21] 




Nts = 3 [0,21,42] 
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Table 5: Comparison of effective masses in the pion (upper half) and rho (lower half) channels, 
averaged over A'ts different source positions which are indicated in the square brackets. denotes 
the performed number of hits per source position. The total number of inversions, A^tot = Nts-N^, is 
constant in each row of the table. 



The results of such an analysis are shown in Table [5| Each row represents a particular 
fixed value of A'tot- In the pion channel there are clear indications that the statistical error 
decreases when more timeslices are used. In the vector channel this effect is less pronounced 
but becomes evident when Nts-N^ ^ 12. Thus, a good strategy to enhance that statistical 
accuracy of mesonic two-point correlation functions is to use as many timeslices as one can 
afford for a fixed total number of inversions. 
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5 Summary and conclusions 



In this paper we have presented a systematic study of random noise sources for the calcula- 
tion of mesonic two-point correlation functions, using 0{a) improved Wilson fermions as our 
discretization. Specifically, we have investigated the effectiveness of the generalized one-end 
trick in the pseudoscalar and vector channels, for pion masses ranging from 450 to 760 MeV. 
Our spatial volumes correspond to a box length of L = 2 fm. The total number of inversions 
of the lattice Dirac operator serves as our measure for the computational cost. 
Our findings are best summarized by a list of empirical observations: 

• At equal computational cost, stochastic volume sources significantly enhance the sta- 
tistical accuracy of correlation functions in the pseudoscalar channel. Here the signal is 
dominated by fluctuations in the gauge configurations, and the extra noise introduced 
by the stochastic procedure is rapidly suppressed. 

• Stochastic volume sources help to stabilize mass estimates in the vector channel. In 
addition to reducing the variance, stochastic volume sources also produce a more pro- 
nounced plateau in the effective mass plot. Compared to the pion channel, however, 
a larger number of hits must be performed, before the improvement is clearly visible. 
In order to observe a clear advantage over ordinary point sources in our studied mass 
range, the numerical effort must at least be doubled for the generalized one-end trick. 

• For fixed numerical cost, the quality of mesonic two-point correlation functions com- 
puted using random volume sources can be improved by averaging over several source 
positions, A'ts > 1, and adjusting the number of hits, A^^ such as to keep the total num- 
ber of inversions fixed. In other words, averaging over more source positions is more 
effective in reducing the variance than increasing the number of hits. 

The use of random volume sources is an attractive method, designed to extract more 
information on hadronic properties from the full propagator matrix at reasonable cost. It is 
particularly useful if the number of available gauge configurations is relatively small. Given 
the convincing performance reported here for the simple case of mesonic two-point functions, 
we have employed random noise sources in our projects to determine the electromagnetic pion 



form factor and the form factors for Ki^ decays 28 
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A Effective mass plots: pseudoscalar channel 
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Figure 5: Effective pion mass plots for the ensembles E2 — (from top left to bottom right). 
Illustrated are the results for the point source (blue squares) and the stochastic volume source (filled 
red cirlces) at fixed cost, i.e. for A^^ = 3 hits. 
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B Effective mass plots: vector channel 
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Figure 6: Effective rho meson mass plots for the ensembles E2 — E5 (from top left to bottom right). 
Illustrated are the results for the point source (blue squares) and the stochastic volume source (filled 
red cirlces) at fixed cost, i.e. for Nr = 3 hits. 
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